Effective Temperatures in Athermal Systems Sheared at Fixed Normal Load 
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We perform molecular dynamics simulations of repulsive athermal systems sheared at fixed normal 
load to study the effective temperature Tl defined from time-dependent fluctuation-dissipation 
relations for density. We show that these systems possess two distinct regimes as a function of 
the ratio T$/V of the granular temperature to the potential energy per particle. At small T$/V, 
these systems are pressure-controlled and Tl is set by the normal load. In contrast, they behave 
as quasi-equilibrium systems with Tl ~ T$ that increases with shear rate at large Ts/V. These 
results point out several problems with using Tl in thermodynamic descriptions of slowly sheared 
athermal systems. 
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Athermal and glassy systems driven by shear or other 
external forces are by definition not in thermal equi- 
librium. Fluctuations in these systems are induced by 
the driving force and do not arise from random ther- 
mal motion. However, there are obvious similarities be- 
tween equilibrium thermal systems and driven systems 
that reach a nonequilibrium steady-state. In fact, there 
have been many attempts to develop thermodynamic and 
statistical descriptions Q of these systems including ki- 
netic theories for driven granular gases 0,the Edward's 
entropy formalism for granular packings 0, Q , and ap- 
plications of equilibrium linear response and fluctuation- 
dissipation (FD)relations |j| to define effective temper- 
atures i n aging 0, Q and sheared glasses 0, @ and com- 
pacting and sheared granular materials [13. H^ . 

The approach that employs FD relations to define ef- 
fective temperatures in athermal and glassy systems has 
shown great promise. First, effective temperatures from 
FD relations for several quantities such as density and 
pressure have been shown to be the same 0, • Sec- 
ond, FD relations can be applied to dense systems with 
elastic particles in contrast to other approaches. FD re- 
lations can also be measured experimentally in athermal 
systems such as foams and granular materials. However, 
it is still not clear whether effective temperatures from 
FD relations can be used in thermodynamic descriptions 
of dense shear flows. Many important questions remain 
unanswered, for example, what variables should be used 
to construct an equation of state and do effective temper- 
ature gradients generate heat flow? We conduct molecu- 
lar dynamics (MD) simulations of sheared, athermal sys- 
tems to begin to address these questions. 

Nearly all numerical simulations that have measured 
FD relations in sheared athermal and glassy systems have 
been performed at constant volume. We will conduct our 
simulations at constant normal load (or external pres- 
sure), instead, and this has several advantages. First, it 
has been shown that sheared, athermal systems behave 



differently in these two ensembles. At constant normal 
load, these s yste ms can dilate in response to an applied 
shear stress |l5l . We will determine how the ability 
to dilate affects FD relations and properties of effective 
temperatures defined from them. Simulations at constant 
normal load will also enable us to determine whether the 
effective temperature is more sensitive to changes in pres- 
sure or shear stress. 

Our principal result is that the effective temperature 
Tl defined from FD relations for density is controlled 
by pressure in slowly-sheared athermal systems. We find 
that for all shear rates below j c , the effective temper- 
ature is independent of shear rate at fixed normal load. 
Above 7c, Tl increases with shear rate and as 7 increases 
further the systems obey quasi-equilibrium FD relations 
at all times. The characteristic shear rate that separates 
the pressure-controlled from the quasi-equilibrium regime 
can be estimated using the ratio of kinetic to potential 
energy [l3|. When the ratio is small, a strong force net- 
work exists and the effective temperature is set by the 
pressure. When kinetic energy dominates, velocity fluc- 
tuations are large, which causes frequent rearrangements, 
diffusion, and quasi-equilibrium behavior. 

We now provide the essential details of the MD sim- 
ulations. The systems contained N/2 large and N/2 
small particles with equal mass m and diameter ratio 1.4 
to prevent shear induced ordering and segregation |17| . 
The particles interacted via one of the following pairwise, 
purely repulsive potentials: 
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where e is the characteristic energy scale and <7y = 
(<7j + <7j) /2 and are the average diameter and separa- 
tion of particles i and j. The repulsive linear (a = 2) and 
Hertzian (a = 5/2) spring potentials (Eq. ^) have been 
used to model granular materials [l8L Eg ]. In contrast 
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to Eq. m the repulsive Lennard- Jones (RLJ) potential 
(Eq. |2l has an infinite repulsive core at small . Both 
interaction potentials are zero for rij > (Tij . 

For athermal dynamics, the position and velocity of 
each particle in the bulk is obtained by solving 



(3) 



where F[ — — ^2jdV(rij)/drijfij is the total repulsive 
force on particle i, Vi is the velocity of particle i, b > is 
the damping constant, and sums over j only include par- 
ticles overlapping particle i. Shear is imposed by moving 
a rough and disordered top boundary in the ir-direction 
at fixed speed u, while a similar bottom wall a distance 
L y away remains stationary. 

To sustain a fixed normal load P ext in the y-direction, 
the top boundary was moved rigidly according to 

r d 2 L„ „„„ „ , ( dL„ \ 
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where M is the mass of the top wall, F™ is the total 
repulsive force acting on the top wall from interactions 
with particles in the bulk, and F ext = P ex tL x L z (P ext L x ) 
is the external normal force applied to the top wall in 3d 
(2d). The equations of motion were solved using a 5th 
order Gear predictor-corrector integration scheme |20j . 
After an initial transient at each u and P ex t, forces in 
the y-direction balance on average, the system fluctu- 
ates about (L y ), and a linear shear flow is established 
with shear rate 7 = u/(L y ). Relatively small system 
sizes N = 1024 were studied to inhibit the formation of 
nonlinear velocity profiles |2l| . Periodic boundary condi- 
tions were implemented in the x (z) direction. We use a, 
e, and er/ym/e, where a is the small particle diameter, 
as the units of length, energy, and time, respectively. All 
quantities are expressed in reduced units below. 

In equilibrium systems, the fluctuation-dissipation the- 
orem requires that the autocorrelation function of a phys- 
ical quantity is proportional at all times t to that quan- 
tity's response to a small conjugate perturbation, and 
the proportionality constant gives the temperature. Re- 
sponse and correlation are not proportional at all times 
in sheared glassy and athermal systems 0. However, 
several recent studies of FD relations for density in these 
systems have shown that they can still be used to define 
an effective temperature that characterizes shear-induced 
fluctuations at long-time scales 0, . Similar measure- 
ments in 3d sheared athermal systems, but at fixed nor- 
mal load, are displayed in Fig. ^ (a) . We plot the inte- 
grated response R(t) of the density of the large particles 
at wavevector k = kr,z 



p(k,t) 



N, 
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to a spatially modulated force with period 2ir/ko applied 
along k at t — and all subsequent times versus the 
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FIG. 1: (a) R(t) vs. C(t) (solid lines) for density p(k) at 
k = 9z for 3d systems with RLJ interactions and b = 0.5 
sheared at fixed normal load P EXt = 0.1. The dotted (dashed) 
lines have slope —1/Tl (— 1/TJsr). Five shear rates are shown: 
7 = 0.53, 0.30, 0.07, 7.7 x 10~ 3 , and 2.0 x 10~ 4 from bottom 
to top. Tl is constant for 7 < 7 C ~ 0.2. (b) Volume fraction 
4> vs. 7 for the same systems in (a). 4> rcp ~ 0.648 for 3d 
bidisperse systems. 

density autocorrelation function C(t) [22]. This figure 
shows that response versus correlation for p(k) can be 
nonlinear at short times in sheared athermal systems. In 
contrast, the slope of R(C) is constant at long times, 
which allows the long-time effective temperature Tl to 
be defined as 



-1/T L 



dR 
dC 



(6) 



Previous studies have shown that Tl does not depend 
strongly on the magnitude @ or direction of k 0. The 
slope of R(C) for p(k) at k = koz in the t — ► limit is 
— I/T5, where T5 = m(v 2 ) is the granular temperature. 

In Fig. (a), we show R(C) for several shear rates at 
fixed normal load P ex t =0.1. At low shear rates, the 
long-time slope of R(C) is independent of 7. However, 
when 7 increases above a characteristic shear rate 7 C , the 
long-time slope of R(C) begins to decrease with increas- 
ing 7. Thus, at constant normal load, Tl is constant for 
7 < 7c but increases for 7 > j c . Fig. [2(b) shows the 
variation of the volume fraction <p with 7 for the systems 
in (a). To compensate for increases in normal force from 
increases in 7, the systems expand in the y-direction. 
Fig. (b) emphasizes that effective temperatures depend 
on the time scale over which they are measured even in 
systems that are 20% below random close packing. 
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FIG. 2: (a) T s and (b) T L vs. 7 for 3d systems with RLJ 
interactions and b — 0.5 sheared at fixed P ext = 0.02 (down- 
ward triangles) or fixed volume fraction (f> — 0.68 (leftward 
triangles). The dotted line in (a) has slope 1.4. The inset to 
(b) shows the internal pressure P vs. 7 at fixed <j> = 0.68. 

The results for R{C) for p(k) at constant normal load 
presented in Fig. ^a) are intriguing; they suggest that 
Tl is fixed by the pressure in slowly sheared systems. 
To investigate this question further, we compare Tg and 
Tl in systems sheared at constant normal load (CN) and 
constant volume fraction (CV) in Fig. |2 We find that 
Ts is insensitive to the choice of the ensemble; Ts scales 
as a power-law in 7 over three decades in both ensem- 
bles. However, the 7 dependence of Tl at CN and CV 
differs significantly. Tl in the two ensembles were ini- 
tially matched at 7 w 10 -4 and P ex t — 0.02 by setting 
(j) in the CV simulation equal to (0) from the simulation 
at fixed P ex t = 0.02. At fixed <fi, Tl increases by more 
than a factor of 3 over the range of 7 studied. The inset 
to Fig.[2](b) suggests that the steady rise in Tl is due to 
the shear-induced increase in the internal pressure P at 
CV. In contrast, Tl at fixed P ex t = 0.02 remains constant 
over the same range of 7. Tl only begins to increase when 
7 > 7c, which in this case is at least 3 decades above the 
quasi-static 7 regime at CV. 

To test the generality of these results, we performed 
measurements of R(C) for p(k) at two values of the 
damping constant (b — 0.5 and 0.04), over a range of 
normal loads from P ext — 0.02 to 10, and over 4 decades 
in shear rate [23j . The results are summarized in Fig. [3] 
Panel (a) shows that for each P ext , there is a wide range of 
7 over which Tl = T® is constant. Tl begins to increase 
above a characteristic shear rate 7 C , which decreases with 
normal load. Fig. [3| (b) shows that Tl begins to deviate 
from T® as Tg increases. In this regime, large velocity 
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FIG. 3: T L plotted vs. (a) 7 and (b) Ts. Panel (c) shows 
Tl/T® vs. the ratio of Ts to the total potential energy per 
particle V. The results were obtained from 3d systems with 
RLJ interactions at fixed P ex t- Five loads P ext = 0.02 (down- 
ward triangles), 0.05 (upward triangles), 0.1 (diamonds), 0.5 
(squares), and 1 (circles) and three loads P ext = 2 (rightward 
triangles), 5 (leftward triangles), and 10 (small circles) were 
studied at b = 0.5 (open symbols) and 0.04 (filled symbols), 
respectively. In panel (a), the dotted line defines the low 
shear rate value T° of the long-time effective temperature at 
Pext = 0.5. In panel (b), Tl = Ts is indicated by the dotted 
line. In panel (c), the dotted line marks the ratio of kinetic 
to potential energy above which Tl /T® begins to increase for 
these systems. 

fluctuations give rise to frequent rearrangement events 
and particle diffusion. At each P ex t, there is a Ts (or 7) 
regime where sheared athermal systems behave as quasi- 
equilibrium systems with response proportional to cor- 
relation at all times and a single effective temperature 
Ts w Tl that increases with shear rate. This nontrivial 
result is not found in sheared glasses because they are 
thermostatted with Ts below the glass transition tem- 
perature and remain arbitrarily far from equilibrium as 
the shear rate is tuned. 

We seek a general criterion to determine whether 
sheared athermal systems exist in the pressure-controlled 
or quasi-equilibrium regime. As a first step, we consider 
the ratio (3 of the granular temperature Ts to the total 
potential energy per particle V . For each value of b, we 
find that there is a range of normal loads P ex t > P^ xt 
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FIG. 4: T° vs. Pext for several sheared athermal systems 
at b = 0.5: 3d RLJ (triangles), 2d RLJ (pluses), 3d linear 
spring (circles), and 3d Hertzian spring (squares). The dashed 
line has slope 0.85. P ext > 10 _1 (10" 2 ) were not considered 
for linear (Hertzian) springs because they caused unphysical 
overlaps. 

where Tl/T® versus Ts/V is independent of P ex t, with 
P£ xt increasing inversely with b. Fig. (c) shows that 
when P ex t > Pe Xt , there is a characteristic ratio f3 c ~ 0.1 
independent of damping that separates the pressure- 
controlled and quasi-equilibrium regimes. These results 
also imply that Tl in sheared glasses is controlled by 
pressure since they have Ts/V <C 1. The behavior at 
large Ts/V in the low-pressure regime P ex t < Pext 1S 
more complex and will be discussed elsewhere pEj . 

In Fig. 0] we plot the long-time effective temperature 
T° in the pressure-controlled regime vs. P ex t for sheared 
athermal systems in 2d and 3d with repulsive linear and 
Hertzian spring and Lennard- Jones interactions. All of 



the data collapse onto a power-law with exponent 0.85 ± 
0.02. If both the normal load [3| and T£ scaled linearly 
with the yield stress of the material, the exponent in 
Fig. 0] would be 1. Our preliminary results indicate that 
the yield stress and pressure are proportional 25] , which 
implies that T° scales sublinearly with the yield stress. 

We used MD simulations to study properties of the 
effective temperature Tl defined from FD relations for 
density in athermal systems sheared at fixed normal load. 
These systems possess two distinct regimes as a func- 
tion of the ratio Ts/V of kinetic to potential energy. At 
small ratios, these systems are pressure-controlled and 
Tl is set by the normal load. At large Ts/V, they be- 
have as quasi-equilibrium systems with Tl ~ Ts that 
increases with shear rate. These results point out several 
difficulties in using Tl in thermodynamic descriptions of 
slowly sheared athermal systems. First, the variables 
Tl, pressure, and density alone do not provide a com- 
plete description of these systems since Tl and pressure 
can remain constant while density can vary substantially. 
Also, these results suggest that when two slowly sheared, 
athermal systems are placed in contact but maintained at 
different pressures, Tl in the two systems will not equi- 
librate. Thus, in this regime Tl does not behave as a 
true temperature variable. We are currently attempting 
to identify a set of variables that can be used in a ther- 
modynamic description of dense granular shear flows. 
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